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(— I We consider a single-server queue with renewal arrivals and i.i.d. 

service times, in which the server employs either the preemptive Short- 
est Remaining Processing Time (SRPT) policy, or its non-preemptive 
variant, Shortest Job First (SJF). We show that for given stochas- 
tic primitives (initial condition, arrival and service processes), the 
model has the same fluid limit under either policy. In particular, we 
conclude that the well-known queue length optimality of preemptive 
SRPT is also achieved, asymptotically on fluid scale, by the simpler-to- 
implement SJF policy. We also conclude that on fluid scale, SJF and 

^\J ' SRPT achieve the same performance with respect to response times of 

► | \ the longest-waiting jobs in the system. 
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1 Introduction 

Schrage's [10J classic result asserts that for given stochastic primitives (initial 
condition, arrival and service processes), the preemptive Shortest Remaining 
Processing Time (SRPT) policy minimizes queue length at all times, over all 
(work-conserving) policies. Recently, Gromoll, Kruk, and Puha [6] showed 
that, in the heavy-traffic regime, the diffusion scaled queue length process 
under SRPT also achieves the optimal lower bound W(t)/x* for all t > 0, 
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where W(-) is the diffusion limit of the workload process and x* is the 
supremum of the support of the service time distribution. 

Under SRPT, preemptive priority is given to the job that can be com- 
pleted first. More precisely, consider a single server queue with renewal 
arrivals and i.i.d. service times, and let I(t) index, in arrival order, those 
jobs that are in the queue at time t. For a job i £ X(t), let Vi(t) denote 
its residual service time at time t, or the remaining amount of processing 
time required to complete this job. If j € X(t) is the smallest index such 
that Vj(t) < Vi(t) for all i € T(t), then under SRPT, £vj(t+) = -1 and 
jj.Vi{t+) = for all i G l(t) \ j. 

Due to its simple description and attractive optimality properties (see 
also Smith [32]), the SRPT policy and its relatives have received much at- 
tention over the past several decades; see Schreiber [TT] for a survey, as well 
as [21 H] and references therein for more recent work. 

Despite its theoretical advantages, SRPT may not always be the best 
choice in practice. Practical disadvantages include inexact job size infor- 
mation, the need to continuously monitor the state of the queue and keep 
this state as a continuum of priority classes, inefficiencies associated with 
frequent switchovers, and long response times for large jobs. See for example 

mm- 

These disadvantages have led researchers to consider a number of vari- 
ants of SRPT, all of which retain the essential feature of prioritizing smaller 
jobs in some way. Care is needed when reviewing the literature, as differ- 
ing terminology is in use. Usually, "SRPT" refers to the preemptive policy 
defined above. This policy may also be referred to as Shortest Time to 
Completion First, or STCF. 

Shortest Job First (SJF) is the non-preemptive variant of SRPT. An 
SJF server still selects the smallest job in the queue. But once begun, 
this job is served to completion before the server makes another selection. 
This policy may also be referred to as Shortest Job Next (SJN). There are 
also preemptive versions of SJF and SJN, which differ from SRPT because 
they prioritize based on initial rather than remaining job size. Nuyens and 
Wierman [13] have introduced the class of e-SMART policies, which include 
the preemptive policies mentioned above as well as policies that allow for 
inexact job size information and finite sets of priority classes based on job 
size. 

In this paper, we compare the preemptive SRPT policy to its non- 
preemptive SJF variant. Our analysis implies that, in many respects, these 
two policies achieve the same performance (although not in all respects; see 



below). Although SJF does not resolve all of the criticisms that have been 
raised for SRPT, it does have important advantages from an implementa- 
tion stand point. In particular, there is no need to continuously monitor the 
state of the queue and it avoids any penalties for interrupting service. Con- 
sequently, practitioners interested in performance measures for which there 
is no essential difference between the two policies may opt for the simpler 
SJF policy, and this paper provides the theoretical justification for such a 
choice. 

1.1 Main result 

Consider an r-indexed sequence of single server queues, each with renewal 
arrivals, generally distributed i.i.d. service times, and random initial state. 
In the rth model, let -ZI(-) denote the measure- valued state descriptor when 
the server employs preemptive SRPT, and let Z r (-) denote the state de- 
scriptor under SJF. These processes track the residual service times of all 
jobs, as they evolve according to each respective policy; see Section [2] for 
details. Let Zp(-) and Z r (-) denote these processes under fluid scaling. We 
will impose the standard asymptotic assumptions (f28]) - ([32|) of Section |U In 
particular, there is a limiting random initial measure Zq, a limiting arrival 
rate a, and a limiting service time distribution v, such that the limiting 
traffic intensity is at most one. 

In this setting, Down, Gromoll, and Puha [4] showed that asr-> oo, the 
sequence {Zp(-)} converges in distribution to a process Z*{-) that is almost 
surely an SRPT fluid model solution for the data (a, v) and initial condition 
Z*(0), where Z*(0) = Zq in distribution; see Theorem 3.9 in [4]. 

An SRPT fluid model solution for (a, u) with initial condition £ is a 
deterministic measure- valued function £(•) satisfying ("(0) = £ and a cer- 
tain family of dynamic inequalities; see Section 2.2, Definition 2.2 in [4]. 
Properties of SRPT fluid model solutions were extensively studied in |3]. 

The main result of this paper is that the sequence of state descriptors 
under the SJF policy has the same fluid limit: 

Theorem 1. Under the asymptotic assumptions (|28p -( l32j) . the sequence 
{Z r (-)} converges in distribution to the measure valued process Z*{-), where 
Z*{-) is almost surely an SRPT fluid model solution for the data {a,v) and 
initial condition Z*(0), and Z*(0) is equal in distribution to Zq. 



1.2 Implications for comparing performance 

Since the queue length processes under SJF and SRPT are given by the total 
mass of the state descriptors Z r {-) and -2I(-) respectively, an immediate 
consequence of the above theorem is that, asymptotically on fluid scale, 
SJF enjoys the same queue length optimality as SRPT. 

A second consequence relates to state-dependent response times. A job's 
state-dependent response time is the time until it exits the system, condi- 
tional on the state of the system (the configuration of all current residual 
service times) when it arrives. The fluid limit Z*{-) can be used to calculate 
a fluid approximation s(x) to the state-dependent response time of a job of 
size x, as a function of an initial state £; see dj. Thus Theorem [1] im- 
plies that, asymptotically on fluid scale, a job of size x arriving to a queue 
in state £ will have the same state-dependent response time under SJF or 
SRPT. Moreover, the system will be in the same state £ under either SJF 
or SRPT. 

Note that SJF and SRPT exhibit differing performance for the mean 
response time of a job of size x (that is, averaged over all possible states 
encountered on arrival when the system is in steady state). Bansal and 
Harchol-Balter [2] show for heavy-tailed service time distributions, that the 
mean response time under SRPT is orders of magnitude lower than under 
SJF, while the very largest jobs have smaller mean response times under SJF. 
See also Harchol-Balter, Sigman, and Wierman [7\ for explicit formulae for 
mean response times under SRPT, SJF, and several other protocols. 

The reason SJF and SRPT have the same fluid approximation for state- 
dependent response times, is that fluid scaling compresses time such that 
individual jobs - even large ones - exit almost instantaneously once they 
begin service. In fact, any arriving job with service time smaller than the 
largest job to begin service during the current busy period, or frontier, will 
have negligible response time on fluid scale. Jobs with service times larger 
than the frontier wait a positive time on fluid scale to begin service for the 
first time. 

Thus, for either policy, the fluid limit shows that response times for jobs 
larger than the frontier are orders of magnitude longer than for those below 
it. In particular, the fluid limit captures the profile of reponse times only 
for those jobs that experience the very longest response times in the system, 
namely, those that are bigger than the frontier when they arrive. And given 
the initial state, these response times are the same under SJF or SRPT. This 
conclusion is consistent with Nuyens and Zwart [8j, where it is shown that 
the steady-state response time has the same asymptotic decay rate under 



SJF and SRPT. 

1.3 Notation 

For real numbers a and b, let af\b = min{a, b} and aVb = max{a, b}. Denote 
the natural numbers 1,2,... by N and the non-negative real numbers [0, oo) 
by R+. The measure that puts one unit of mass at x £ R+ is written 8 X , 
and the meaure 5+ is 5 X if x > and the zero measure (denoted 0) if x = 0. 
Let M denote the Polish space of finite, nonnegative Borel measures on 
R+, endowed with the weak topology [9J. A metric inducing this topology 
is the following Prohorov metric on M. For A C R+, let A e = {y G R + : 
inf xGj 4 \x — y\ < e}. Then for £, £ € M, define 

d[£, C] = mf{e > : £(B) < ((B e ) + e and C(S) < £(B e ) + e 

for all closed B CR+}. 

For C € M and a £-integrable function g : R + — >• R, define (5, £} = 

J R g(x)C(dx). Let x( x ) = ^ be the identity function, and for any set A, 

let 1,4 denote its indicator function. For a real- valued function /, write f + 

and f~ for its positive and negative parts respectively. Unless otherwise 

specified, all processes are assumed to be right continuous with finite left 

limits (rcll). In particular, all measure- valued processes in this paper take 

values in the Skorohod space D([0, 00), M) endowed with the Ji-topology. 

p 
Convergence in probability is denoted — >, weak convergence of elements 

of M is denoted — ►, and convergence in distribution of random objects is 

denoted =>. We adopt the convention that a sum of the form YULn w ith 

n > m, or a sum over an empty set of indices equals zero. Finally, we also 

define min = inf = 00. 

2 The model: SJF and SRPT 

This section lays out the precise definitions of our model, including the two 
policies under consideration. We will compare the SJF and SRPT policies 
pathwise, using the same sample space (0, J 7 , P) and same stochastic prim- 
itives. The stochastic primitives consist of an exogenous arrival process, a 
sequence of service times, and an initial condition. 

Stochastic primitives. The arrival process E(-) is a rate a £ (0, 00) 
(possibly delayed) renewal process. For t > 0, E(t) is the number of jobs 
that have arrived to the queue during (0, t]. The initial service times of these 



jobs are taken from an i.i.d. sequence {vjj : j = 1, 2, . . . } of strictly positive 
random variables, with distribution v having finite mean. Define the traffic 
intensity p = a(x, v). 

The initial condition is given by a non-negative integer-valued random 
variable Zq with finite mean, and a sequence of strictly positive random 
variables {wj : j = ...,— 2, — 1}. The random variable Zq represents the 
number of jobs in the queue at time 0, and their initial service times are 

{w- Zo ,---,w-i}- 

Define the random index set J = {j 6 Z : j > —Zq} \ {0} and, for each 
t > 0, define the random index set «/(£) = J n {—Zq, . . . ,E(t)}. Then J 
indexes all jobs that will ever be in the system, and </(£) indexes all jobs 
that have been in the system by time t. Note that neither set ever contains 
the index 0, and that J(t) = if Zq = E(t) = 0. Note that the present 
indexing differs slightly from that used in [3], but this does not affect our 
analysis. 

Workload. Define a measure- valued load process V(-) and a load pro- 
cess V(-) by 

E(t) 

V(t) = ( x ,V(t)) , 

for all t > 0. Define the workload process W(-) and cumulative idle time 
process /(•) by 



W(0)= £ 
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I(t) = sup[W(0) + V(s)-s}-, y ' 

s<t 

W(t) = W(0) + V(t) - t + I(t). 

Next, we define the SJF and SRPT policies, as well as the corresponding 
state descriptors that describe the evolution of the system under each policy. 

Shortest Job First. In order to define the residual service times under 
SJF, we first define the succesive times {7™} at which the server begins 
processing a new job. Let 70 = ?/o = jo = and define the set Jq = 0. Now 



define inductively, for n = 1, 2, . . . , 

7„ = inf{t > 7„_i + y n _i : Z + E(t) > n}, 

y n = mm{wj : j G J(j n ) \ J n -i}, , 

j n = min{j G J(7„) \ J n _i : w,- = y n }, 

•Jn — "n— 1 U iJrtj"' 

For each n = 1,2,..., the time 7„ is the nth time at which the server begins 
processing a new job, y n is the service requirement of this job, j n is its index, 
and J n is the set indexing all jobs that have begun service by time j n . For 
each j G J, let 

n(j) = min{n € N : j G J ra }. 

Then n(j) is the overall order in which job j is served, and the start time of 
job j may be defined as 7 n (j) (with 700 interpreted as 00, which may occur 
for example in a supercritical model in which job j never begins service). In 
particular, note that y n (j) = w j f° r an j € J such that 7 n (j) < °°- At time 
t > 0, define the attained service of job j € J as 

Sj (t) = (t_ 7n(j) )+ Au)j, 
the cumulative service provided as 

5(t)= 2 S M 

jeJ(t) 

and define the residual service time of job j G J as 

Wj(i) = Wj — Sj(t). 

We now define the state descriptor under SJF as the measure-valued 
process 

Finally, for each t > 0, let -F(t) denote the largest job started by the SJF 
server by time t. That is, 

F(t) = max {y n : 7„ < t} . 

neN 

Shortest Remaining Processing Time. For job j G J, let fj(t) 
denote its residual service time at time t when the server employs the SRPT 



policy; we refer the reader to [1] for a detailed definition. Then the measure- 
valued state descriptor under SRPT is given by 

j'6J(t) 

Throughout the paper, the subscript p for "preemptive" is used to distin- 
guish performance processes under SRPT from their analogues under SJF. 
In [3], the left edge process of the state descriptor Z p {-) is defined as 

L p (t) = sup{x G R+ : (1^.) ,Z p (t))=0}, t > 0. 

The current residual service time process is then defined, for all t > 0, by 

I 0, otherwise, 

and the SRPT frontier process is defined by 

F p (t) = sup C p (s), t > 0. 

0<s<t 

3 Pathwise comparison 

In this section, we establish several results that relate SJF to SRPT on a 
pathwise basis. The main goal is to show that SJF and SRPT have the 
same frontier process almost surely. This will provide the key element in 
the proof of Theorem [1] in Section [U 

First, we show that the sum of the residual service times under SJF, as 
defined above, equals the workload defined in (P). 

Lemma 1. Almost surely, for all t > 0, S(t) = t — I(t), and W(t) = 

Proof. Clearly S(t) < W(Q) + V(t) for all t > 0. Recall that a regular point 
of a function is a point at which it is differentiable. We first show that, for 
all regular points t, 

S(t) = if and only if S(t) = W(0) + V{t). (3) 

Note that almost surely for all j G J, the derivative Sj(t) G {0,1} for all 
regular points t of Sj(-), and Sj(t) = 1 if and only if t G (in(j)i^n(j) + Un(j))- 



Further, the definitions above imply that the finite elements of {n(J) : j G J} 
are distinct and the time intervals {(7 n ,7n + Vn) '■ n G N} are all disjoint. 
This implies that S(t) G {0, 1} for all regular points £. (Note that if £ is a 
regular point of £*(•) but not of all Sj(-), j G </(£), then we must have S(t) = 1 
anyway). Moreover, S(t) = if and only if t < 71 or £ G ( r y n + y n , 7 n +i) for 
some n G N. If £ < 71, then S(t) = = W(0) + V(t), and the converse is also 
true. If £ G (7 n + y n , 7 n +i) for some n G N, then (J2j) implies Zo + E(t) = n, 
which implies that Sj(t) = Wj for all j G J(£). Thus, S(t) = W(0) + V(t), 
and the converse is also true. This proves Q. 

Observe that almost surely, S(-) is continuous, starts at zero, has deriva- 
tive in {0,1} at regular points, is bounded above by W(0) + V(-), and 
satisfies ([3]) at all regular points t. Since W(0) + V(-) is right-continuous 
and piecewise constant, this uniquely determines S(-) almost surely. 

As is well known, I(t) G {0, 1} at all regular points t, and I(t) = 1 if 
and only if I(t) = (W(0) + V(t) - t)~ and 1(f) > 0. So f t (t - I(t)) = if 
and only if £ - 1(f) = W(0) + F(t). Also, t - I(t) < W(0) + V(t) for all 
£ > 0. Since • — /(•) is continuous and starts from zero, it is also uniquely 
determined by W(0) + V(-) almost surely. We conclude that S(t) = £ — I(t) 
for all £ > almost surely. 

The second part follows directly from the first part and equation (jTJ), 
since 

W(t)=W(0) + V(t)-t + I(t)= Yl Wj-S(t)= Yl Wi(t). 

ieJ(t) jeJ(t) 

D 
Note that by Lemma [H for all £ > 0, 

W(t) = ( X ,Z(t)). (4) 

Next, we prove a result describing the workload at or above the frontier 
under SJF, at certain random times. To that end, define the random set 

J = {j G J : 7 n(i) < 00 and Wj = F(j n ^)}. 

These are frontier jobs, with service time equal to the frontier when they 
are started. The set of start times of frontier jobs is T = {infj) '■ 3 £ J}- 
Let {tj, : k G N H [0, \J~\]} be the ordering of T such that r^ < r^ +1 for all 
k < \J"\. Note that J is an infinite set in most settings of interest. But it 
can happen that \J\ < 00, for example when p > 1 or when there are initial 
jobs larger than the supremum of the support of v. 



Lemma 2. (i) Almost surely, for all k G N with k < |J7"|, 

( X l[F(T h ),oo),Z(T k )) = W(T k ). (5) 

(ii) Almost surely, for all k < \<J\ and t G (Tk,Tk+i], 

E(t) 
(xMF(r k ),oo),Z(t)) =W(T k )-F(T k )+ Yl W i 1 {w i >F{T k )}- (6) 

j=E(r k )+l 

Moreover, if k = \<J\ < oo and t G (T k ,oo), then © holds as well. 

Proof, (i) Fix k G Nfl [0, \J~\]. By definition, r k = 7 n (j ) f° r exactly one jo G 
J, and F(r k ) = F(7 n y )) = w J0 . Consider any j G J (7}.). If 7n(j) < 7*, then 
([2]) implies that 7 n (j) + Vn(j) < 7n(j"o)> an< ^ so w j( T k) = 0. On the other hand, 
if 7n(i) > Tfe, then u/,-(r fc ) = 10, and (|2j) implies that j G J(ln(j )) \ J n(j )-i- 
In particular, iwj > y n (j ) = Wj . So at time Tfc, all arrived jobs have residual 
service time either equal to zero or bounded below by F(r k ). This implies 
that the first right-hand term in 

W(T k ) = (x\o,F(r k )),Z{T k )) + (xl[F(T h ),co),Z(T k )) 

equals zero and proves ©. 

(ii) To prove ([6]), fix k < \J\ and t G (7^,7^+1]. Then 

{xl[F( Th ),oo),Z(t)) = Y, Wj(t)l {wj{t) > F{Tk)} . (7) 

jeJ(t) 

Observe that Wj (t) < Wj = F(r k ), because t > 7 n (j ) an d ln(j ) ^ T. So 

W Jo( t ) 1 {w J0 (t)>F(r k )} = 0, (8) 

while 

w Jo( T k)l{ Wn ( Tk )>F(T k )} = F(r k ). (9) 

Now consider j G J(t) \ jo- We will show that 

W j( t ) 1 {w J (t)>F(r k )} = w j( T k)l{ Wj (r k )>F(T k )}- (10) 

Case one: If 7 n (j) < T ki then Wj(r k ) = by ([2]). Since Wj(-) is nonin- 
creasing and t > T k ,Wj(t) = as well and so (TlOl) holds. 

Case two: If 7 n y) G (rfc,Tfc + i), then j £ J and so lo, < F(7 n (j))- 
Since F(-) is right-continuous, piecewise constant, and can only jump at 
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(a subset of) times r G T, this implies Wj < F(t}.). Then lf w (t)>F(T k )} = 
li w . (r k )>F(r k )} = an d equation (fTUj) holds. 

Case three: If J n (j) > r fe+i) then "i n u\ > t and so Sj(t) = 0, that is, 
Wj(t) = Wj(Tk) and equation (fTOj) holds. 

This proves (fTUj) . Combining (J7D with (jSJ) and (fTU|) and noting that if 
j > E(jk) then Wjir/.) = Wj, we obtain 

E(t) 

(xl[F(r fc ),oo),^(i)) = J] W i( T fc) 1 Wr fe )>^(T- fc )}+ X] "^{wj^fa)}- 

jeJ(T k )\J0 j=E(T k )+l 

By Q, including the j'o-term in the first sum is compensated by subtracting 
F{jk), which yields 

E{t) 

(x\F(r k ),oc),Z{t)) = (xl[F(T k ),oo),Z(Tk)) ~ F{T k ) + ^ W 3 l {wj>F{r k )}- 

j=E(r k )+l 

Applying ([5]) proves ([6]). 

It remains to show © when k = \J\ < oo and t E (t^, oo). The argument 
is identical except for Cases two and three in the proof of (flOl) . Case two 
becomes 7 n (j) £ (^,00). Then j ^ J and u>j < F{^ n ij-\). Since F(-) is 
constant on [r^oo), Wj < F{jk) and (fTO]) follows as before. Case three 
becomes 7 n (j) = oo. Then Wjfa) = Wj(t) = Wj, which implies (fTUj) , The 
remainder of the proof follows from (jT|)- (|10p as before. □ 

We need a similar result for SRPT, proved in the next two lemmas. For 
each j £ J, the SRPT start time of job j is defined rjj = inf {t > : Vj (t) = 
-1}. Let 

J p = {j G J : rjj < oo and Wj = F p (r)j)}, 

let T p = {rjj : j G J p }, and let {dk '■ k € N n [0, \Jp\]} be the ordering of T v 
such that Gk < Cfc + i for all k < \Jp\. 

Lemma 3. Almost surely, (i) if t G 7~ p , then C p (t) = F p (t), and (ii) if 
C p (t) = F p (t) > 0, then t G T p . 

Proof. Fix t G T p , so t = rjj for some j G J p . Then Wj = F p (r]j). Evidently 
Vj(t) = Wj > 0, so Z p (t) ^ 0, and therefore C p (t) = L p (t). Now suppose 
(l[o )W ), Z,p{t)) > 0. Then there is at least one job in the queue at time t with 
residual service time y G (0, Wj), contradicting the assumption that SRPT 
starts job j at t. This shows that C p (i) > F p (t). Since F p (t) > C p (t) by 
definition, this establishes (i). 
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For (ii), fix t such that C p (t) = F p (t) > 0. By g] Lemma 5.3(i), 
(l[o,c (t)))-2p(0) = almost surely. Additionally, noting differences in the 
indexing of jobs, [1] Lemma 5.3(h) implies that for all j G J such that 
iuj > F p (t), Vj(s) = Wj for all s G [0, t]. Since C p (i) > 0, Z p {t) ^ and 
so (l[C p (t),oo)>-Zp(*)) > 0. Thus the set {j G J(t) : C p {t) = Vj(t) = Wj} is 
nonempty. Let j* be the smallest element of this set. Almost surely, there 
exists e G (0, C p (i)) such that £"(£ + e) = !?(£), that is, no new jobs arrive 
during (t,t + e\. Then £v r (s) = on [0,t) and £vj*(s) = -1 on (t,t + e), 
so r/j. = inf{s > : Vj*(s) = —1} = t. Thus wj* = F p (t) = F p (j]j*) and 
hence j* G Jp and t £ T p - □ 

Lemma 4. (%) Almost surely, for all k G N «;i£/i A; < |j7pj ; 

(xl[F P K),oo),^ P (^fc)> = W(o-fc). (11) 

(wj Almost surely, for all k < \Jp\ and t G (cfej^fe+i]; 

E(t) 

i=E(<r fe )+l 

Moreover, if k = \Jp\ < oo and t G (a*;, oo), £/ten (|12p ZioZds as well. 

Proof (i) Fix k G Nn [0, |J p |]. By @] Lemma 5.3 (i), (l %Cp{t)) ,Z p {t)) = 
for all t almost surely. Since a^ G T p , F p (ak) = C p {ak) by Lemma and it 
follows that (1[o,f (a k ))-,Z p (ak)) = 0. This implies that the first right-hand 
term in 

W(a k ) = (xlio,F p („ k )),Z p (o- k )) + {x\F P { ak ),oo)iZ P {ok)) 

equals zero and proves (fTTI) . 

(ii) To prove f|12|) . fix k < \J P \ and fix t G (a/-, Pk+i]- Then 

{x\F p {a k ),oo),Z p {t))= ^2 Vj{t)l {vj(t) > Fp(<7k)} . (13) 

Let jo G J p be the job started at er^.. Then u;j = F p {ak) and i>j (i) < Wj 
since t > at = i]j . Thus, 

v Jo(t) 1 {v JO (t)>F p (<r k )}=0, (14) 

while 

tyoC *) 1 ^,,^*)^.^)} = Fpivk)- (15) 
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Now consider j G J(t) \jq. We will show that 

V j( t ) l {v 3 {t)>F p {a k )} = V j( CF k)^{v J (a k )>F p (a k )}- (16) 

Case one: If rjj < a k , then Vj(t) = — 1 on {r]j,r}j + e) for some e > 0, so 
v j{?k) < Wj < F p (rjj) < F p (<Jk). Since Vj(-) is nonincreasing and t > a k , we 

have l^t^K)} = ^lyfoj^fa)} = °> and CED holds. 

Case two: If rjj G (cfc,<7fc + i), then j ^ J p and so Wj < F p (rjj). Since 
-Fp(-) is constant on [o"fc,<Tfc+i), this implies i>j(£) < Wj(cfe) < u>j < ^p(cfe)- 
Then l{« 3 -(t)>i? p (<T fc )} = ^K^^Ok)} = and (US]) holds. 

Case three: If 77^ > Ofc+i, then ^ > t and so Vj(i) = Vj(ak) and (|16p 
holds. 

This proves (fi"6|) . Combining (fl"3j) with (fT4l) and (fTBT) and noting that if 
j > E{ak) then Vj(a k ) = Wj, we obtain 

E{t) 
(xl[F p {<T h ),oo),Zp(t))= Yl V j(. a k) 1 {v j (<T k )>F p (<T k )}+ Yl W J 1 {^j>F P (a k )}- 

ieJ(<r fc )\io j=£:(<7 fe )+i 



By (I15p . including the jo-term in the first sum is compensated by subtracting 
F p {ak), which yields 

E(t) 
{x\F p {a k ),oo),Zp{t)) = {xl[ Fp {a k ),oo), Zp{a k ))-F p (a k )+ Y W 3 l {wj>F p {cr k )}- 

j=E(u k )+l 



Applying (|Tlj) proves (IT2|) . 

It remains to show (|12p when A; = \J p \ < 00 and £ G (cfc,oo). The 
argument is identical except for Cases two and three in the proof of (116p . 
Case two becomes rjj G (<7fc, 00). Then j ^ J^ and Wj < F p {rjj). Since 
-Fp(-) is constant on [cr fe , 00), u>j < Fp(a k ) and (fT6|) follows as before. Case 
three becomes r/j = 00. Then Vj(a k ) = Vj(t) = Wj, which implies (fTUj) . The 
remainder of the proof follows from (|13p -(|16 p as before. D 



We now show that SJF and SRPT have the same start times of frontier 
jobs, and coinciding frontiers at those times. This is the principal technical 
result of the paper. 

Lemma 5. Almost surely, T = T p and F(t) = F p (t) for all r G T . 

Proof. Clearly, a\ = inf{t > : Z + E(t) > 1} = 71. For the SJF queue, 
w ji = 2/1 = -^(71)) so 3i G J and thus 71 G T, implying t\ = 71- Conclude 
that T\ = a\. Moreover, F{j\) = xmn{wj : j G J(ti)} = F p (tx). Now fix 
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k < \J\ A \J P \ and assume that 77 = 07 and F(t{) = F p (ti) for all / = 1, . . . , k. 
We will show that r k+1 = a k+1 and F(r k+1 ) = F p (r k+1 ). 

Note first that by definition of the start times, Z(jk+i) 7^ and 
Z p (a k+ i) y^ 0. Since both queues have the same workload process W(-), this 
implies that W (r k+ i) AW (a k+ i) > 0, and so Z(a k+ i) ^ and Z p {r k+ i) ^ 0. 
We first show that r k+ \ = a k+ \. 

Suppose that a k +i > T k+ \. Then r k+ \ G {a k , Cfe+i), so by LemmaSland 
the induction hypothesis, 

E{r k+1 ) 
(xllF p (* k ),oo),Zp(T k +i)) = W{jk) -F{r k )+ *Yj W J 1 {^ J >F(r k )}- 

j=E(r k )+l 

Using t = r k+ i in (0 of Lemma [2J rewrite the right side of the display above 
to obtain 

(xl[F p (ff fc ),«0>-Zp(7fc+l)) = (xl[ir( Tfe ) :00 ),^(r fe+ i)). 

Since F(r k ) < F(r k+1 ), the right side is bounded below by (xl[F(T fc+1 ),oo)> ^(tjh-i)), 
and bounded above by W(r k+ i), both of which are equal by ([5]) of Lemma 
EJ Thus, 

(xV P fo),oo)>^»(7*+l)) = W(n+i), 
which implies that 

(l[0,F p (ff k ))>Zp(Tk+l)) = (xl[0,F p (a k )), Z p(l~k+l)) = 0. (17) 

Since Z p (r k+ i) ^ 0, (fTTl) implies that C p (rfc + i) > F p (a k ). But -F p (-) is 
constant on [a k ,a k+ i), so C p (rfc + i) > F p (rfc+i), which implies by Lemma[3] 
that Tfc+i £ 7p, a contradiction. 

Now suppose that <7fc + i < Tfc+i. We proceed analogously to the above 
argument, with the roles of T k+ i,a k+ i and Lemmas [2] and H] reversed. We 
have a k+ i £ (r k ,T k +i), so by Lemma [2] and the induction hypothesis, 

E((T k+1 ) 

(xl[F(r fc ),oo), Z(a k +l)) = W{a k ) - F p (a k ) + ^ W j 1 {w ] >F p (a k )}- 

j=E(a k ) + l 

Using t = cr k+ i in (I12p of Lemma U rewrite the right side of the display 
above to obtain 

(xl[F(r fc ),<x>)>- 2: (c r fc+l)) = (xl[F P K),oo)) 2 p(c r fc+l)}- 
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Since F p (a k ) < F p (a k+1 ), the right side is bounded below by (x^[F p (<r k+1 ),oo), Z P {°k+i)), 
and bounded above by W(a k+ i), both of which are equal by ([TT]) of Lemma 
H Thus, 

(xl[F(T k ),oo),Z(a k+1 )) = W(a k+1 ), 

which, since F(-) is constant on [r k ,T k+ i), implies that 

(xl[0,F(<T k+1 )),Z(cr k +i)) = 0. (18) 

Since Z(a k +i) 7^ 0, ([15]) implies that the set I = { j G J(<7fc + i) : u)j(<jfc+i) > 
^(crfc+i)} is nonempty. Let j* = argmin{7 ra (j) : j G /} be the job in / with 
the earliest start time. Then by (fT8|) . u>j , .„> _ 1 (fJ k +i) = 0, which implies that 
ln(j*)~i+y n {j*)-i < °"fc+i and Z Q +E(a k+ i) > n(f). Conclude from §2§ that 
JnU*) = °"fc+i- Moreover, the definition of F(-) implies that Wj* = F(a k+ i), 
and so a k+ \ G T, a contradiction. 

The previous two arguments imply that T k+ \ = a k +\. We now show that 
F(r k+ i) = F p (T k+ i). Suppose that F p (r k+ i) > F(r k+ i) so that the interval 
[F (r k+ i) , F p (T k+ i)) is nonempty. By definition of T, there is a job j* G J 
such that 7 n (j„) = T k+ \ and Wj t = F(r k+ i). Then Wj t {T k+ i) = Wj t and so 

(x~L{F(T k+1 ),F p (T k+1 )), z ( T k+i)) > 0- (19) 

Note that for j G J{T k+ \), Wj > F p (T k+ \) implies that Wj > F(r k+ i) 
and so Wj(r k+ i) = Wj > F p (T k+ i). Conversely, Wj(T k+1 ) > Fp(r k+ i) clearly 
implies Wj > Wj(r k+1 ) > F p (r k+1 ). So 

Wj(T k+1 )l {Wj{Tk+l) > Fp{Tk+l)} = Wjl {wj > Fp{Tk+l)} . (20) 

Similarly, Wj > F p {T k+ i) implies Wj > F(r k+ i) > F(r k ) = F p (r k ), and so 
Vj(r k+1 ) = Wj > F p (r k+ i). Conversely, Vj(r k+ x) > F p (r k+1 ) gives Wj > 
Vj(T k+ i) > F p (r k+ i) and so 

Vj{Tk+l)l{v 3 {T k+1 )>F p {T k+1 )} = w j l {w 3 >F p {r k+1 )}- (21) 

Together, ([20]) and ([21]) imply that 

Wj(T k+1 )l {Wj{Tk+l) > Fp{Tk+l)} = Y^ V j( T k+l) l {v,(r k+1 )>F p (r k+1 )}, 



j£j(r k+1 ) iGJ(r fe+ i) 



and so 

(xl[F p (T k+1 ),oo),Z(i~k+l)) = (xl[F p (T k+1 ),oo), Z P (T~k+i)} ■ (22) 
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By mi, 

(xl[F(T k+1 ),oc),Z(.T k +l)) = (x\F{r k+1 ),F p {T k+1 )),Z{T k+1 )) 

+ (X^[F p (T h+1 ),o C ),2(Tk+l)) 
= (xl[F(T h+1 ),F p (T k+1 )),Z(T k+ l)) 

+ (xl[F p (r k+1 ),oo), Zp(T~k+l)} ■ 

But the left side and second term on the right side of (|23p are both equal 
to W(rfc+i) by Lemmas [2] and HI So the first term on the right side equals 
zero, contradicting (|19p . This argument is symmetric if the roles of F(rk+i) 
and F p (Tk+i) are reversed, and we conclude that F(rk+i) = F p (Tk+t)- 

This completes the proof if \J\ = \J P \. But this is the only possibility, 
for if \J\ > \J P \, then for k = \^T P \, there exists t^+i G J such that t^+i € 
(<7fc,oo). This yields a contradiction by the argument in the paragraph 
leading to (117p (substituting oo for o~k+i there). A similar contradiction is 
implied by \J\ < \J P \. D 

Corollary 1. Almost surely, F(-) = F p (-) and for all t > 0, 

(xl[F(t),oo),Z(t)) = (xl[F(t),oo),Zp(t))- 

Proof. The first statement follows from Lemma [5] since both F(-) and F p (-) 
are right-continuous and constant on the intervals [t/^t/c+i) for all k < \J\ 
and [tj., cxo) for k = \,J\. The second statement follows from this fact as well, 
combined with the first statement, and Lemmas [21 01 and[5j □ 

Since the frontier processes are identical, the notation F(t) will be used 
hereafter for the frontier process under either policy. We will use the fact 
that the state descriptors under either policy are identical above the frontier. 

Corollary 2. Almost surely, for all measurable f : M+ — > M+ and all t > 0, 

{f l [F(t),oo): Z { t )) = {f l [F(t),oo), Z p{t))- 

Proof. Fix t > 0. Clearly for all j G J(t), Wj(t) > F(t) if and only if 
Wj > F(t) and Wj(t) = Wj, and Vj(t) > F(t) if and only if Wj > F(t) and 
Vj(t) = Wj. So for all measurable /, 

(fl(F(t),oo),Z(t)) = Y^ f( w j) 1 {w j >F(t)} = (fl(F(t),oc),Zp(t))- (24) 

jeJ(t) 

Using the special case f = X an d combining with Corollary [T] yields 

(l {F(m ,Z(t)) = (l {F(m ,Z p (t)). (25) 

Combining ()24|) and (|25[) completes the proof. D 
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We will also need the following simple bound for the workload below the 
frontier in the SRPT queue. By the workload equation in ([I]), 

W(s) = W(t) - (V(t) - V(s)) +t-s- (I(t) - I(s)) 
<W(t) - (V(t) - V(s)) + t - s, 

almost surely for all t > s > 0, since the cumulative idle time process /(•) is 
nondecr easing. By Q, and by splitting the workload at times s and t into 
residual service times that are strictly below the frontier, and those that are 
at or above the frontier, 

(xl[0,.F(»)>-2p(s)) + (xl[F(s),oo),^pO)) 

E{t) 

< {x\o,F{t)),Zp{t)) + (xl[F(t),oo)>Zp(t)) ~ ^2 w j+t- s - 

j=E(s)+l 

Note that Vj(u) > F(u) implies Vj(u) = Wj for all u > and j € J(u). 
Thus, 

(xl[0,F(s))>^ : p( s )) + 22 W J l {vj(s)>F{s)} 

je.J(s) 

E{t) 

< (xl[o,F(t)),Z P (t)) + J2 w J 1 {vj(t)>*'<t)} ~ ^2 Wj + t-s- 
jeJ(t) j=E(s)+i 

Since F(-) is non-decreasing and residual service times Vj(-) are non- increasing 

w 3 l {v 3 (t)>F(t)} < ^1{^.( S )>F( S )} for all j e J(s). Also, Wjl{ Vj (t)> F (t)} < Wj 
for allj = E(s) + l,...,E(t). So 

E(t) 

J2 w J 1 {v J {t)>F(i)} ~ Yl w i ~ Yl w ^{vj(a)>F(s)} ^ °' 

jeJ(t) j=E(s)+i jej(s) 

and thus almost surely for all t > s > 0, 

{xl[o,F{s)),Zp(s)) < (x i -[o,F(t)), 2 p(t)) +t-s. (26) 

4 Invariance of fluid limit 

Now we shift our point of view from a single model to a sequence of models. 
After defining the sequence and making appropriate asymptotic assump- 
tions, we will use the work of the previous section to show that the models 
have the same fluid limit under SJF and SRPT. 
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Let 1Z be a sequence of positive real numbers increasing to infinity. For 
each r £ 11, there is an associated stochastic model with initial condition 
Zq and {w V j : j = ...,— 2,— 1}, arrival process E r (-) and service times 
{w'j : j = 1, 2, . . . }. These primitives have parameters a r , v r , and p r , and 
are defined on a probability space (il r , J-~ r ,P r ). Then for each r £ 1Z, the 
stochastic primitives give rise to measure- valued state descriptors Z r (-) and 
-23(-) under the SJF and SRPT policies respectively, as well as a common 
frontier process F r (-). Define fluid-scaled versions of the following processes: 



F r (t) - 


--F r (rt), 


V{t) -- 


~- -V'(rt), 
r 


V r (t) - 


- \v r ^l 


Z r (t) - 


~- ~Z r {rt), 
r 


z;(t) - 


- ^(rt), 


W r {t) = 


-- -W r {rt) 
r 



(27) 



Let a > and v € M be a probability measure that does not charge 
the origin, such that p = a(x, v) < 1- We impose the following asymptotic 
assumptions, as r — > oo, on the sequence of models. For the sequence of 
exogenous arrival processes, assume that 

-E r (rt)^a(-), (28) 

r 

where a(t) = at for all t > 0. For the sequence of service time distributions, 
assume that 

v r — > v and {v T : r € 7Z} is uniformly integrable. (29) 



It follows from (|29j) that p r — > p. 

Additionally, assume that as r — > oo, 

(Z r (0),( X ,Z r (0)))^(Z ,( X ,Z )), (30) 

where Zq € M is a random measure satisfying 

(x,Zo) < oo a.s., (31) 

(1 {0} ,Z )=0 a.s. (32) 
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It is shown in [3] that, under the assumptions ([SH ]) — ff32]) asr-> oo, there 
is joint convergence 

(V(-),V r (-),W r (0),Z r p (-)) =* (V*(-),V*(-),W*(0),Z*(-)). (33) 

Here, V*(t) = atv for all t > 0, and Z*{-) is almost surely an SRPT fluid 
model solution for data (a, v) and initial condition Z*(0), equal in distribu- 
tion to Zq\ see Theorem 5.16 in [3] as well as Section 2.2 in [3] for a definition 
of the fluid model solutions. Such fluid model solutions are analyzed in detail 
in glEJ. 

Within the proof of (f33j) . it is assumed by invoking the Skorohod rep- 
resentation theorem that all random elements are defined on a common 
probability space (O*, J 7 * ,~P*) such that, almost surely as r — > oo, 

(V(-), V r (-), W r (0),Z r p (-))^{V*(-),V*(-), W*(Q),Z*(-)) , (34) 

uniformly on compact time intervals, where all objects with a tilde denote 
Skorohod representations; see (84) in [4J. As part of this proof, it is shown 
that almost surely, 

fim (xl [0: ^ (f)) ,-^(t)) = 0, for all t > 0, (35) 

where F r (t) = sup sG r ot i sup{x G M + : (l[o )X ),^(s)) = 0} with sup0 = 0. 
(See (96) in the proof of Lemma 5.23(i) in |3j. Note that since p < 1, 5.23(i) 
is the relevant statement. Also, see the beginning of Section 5.3 in [4J; a 
fixed uj is chosen from an event of probability one in £1*. Since Lemma 5.23 
is proved for this fixed uj, (j35|) holds almost surely.) We will use this fact 
to establish a similar statement for the original process Zp(-), but we need 
to upgrade the pointwise convergence to uniform convergence on compact 
time intervals [0,T]. This is shown in the proof of the following lemma. 

Lemma 6. For all T G [0, oo), 



sup {xhnFr {t)) ,Z r Jt)) — ^0, asr^oo. 

te[o,T] x ' 



Proof. Since the Skorohod representations in (|34p are equal in distribution 
to the original processes, it suffices to show that, P*-almost surely, 



>OG 



te[o,T] 



lim sup (xhoFrM^Zpit)) = 0. (36) 
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To that end, let D C O* be the event of probability one on which (|35p holds, 
and for each r G 1Z define the event 



for all < s < t < oo \ . 



Since the dynamic inequality (|26p is true almost surely, it also holds under 
fluid scaling P r -almost surely for each r G 7£, and therefore P*-almost surely 
for each Z r p {-), r G K. Thus, P* (V n f] r eTl £r ) = L 

Fix T G [0,oo) and e > 0. Let = ao,ai, ...,a& = T be a partition of 
[0,T] such that max^^a.,- — aj-i < e. For £ G [0, T], let aj = minjaj : aj > 
t}. Then on T>, 



and on Hre^ 1 ^' ^ or a ^ r e ^ anc ^ * e [0>^1> 



tna> 



- W^W 1 ^))'^ ^) + e - 



So on P n n re7 ^ r , 

limsup sup (xlr 0j pr m) ,2p(t)) < e, 

r->oo te[0,T] X L ' W; X 

which proves (f36|h D 

Next, we establish a bound for the fluid scaled mass near the origin, 
which is valid under either policy. 

Lemma 7. Let T < 00. For each e, 77 G (0, 1) i/iere exists 5 > suc/i i/iai 



liminf P r sup «l [0)<5] ,.Z r (t)> V (l [0)<5] ,^(t))) < e > 1 - r,. 

Proof. Fix e, 7/ G (0, 1). By Lemma 5.11 in |3], there exists <5i > such that 

liminfP r ( sup (l [0>Sl] ,Z;(t)) < e ) > 1 - J. 
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By ([32]) . there exists 62 > such that 

P((l [o , 52] ,2o)<0>l-|. (37) 

Since v puts no mass at the origin, there exists 83 > such that 

aT(l M3] ,i/)<|. (38) 

Let 5 = min{<5i, 62, £3} and define the events 

n£ = i sup (l M ,4(i))<el, 
[te[o,T] J 

^ = {{1 [0)5] ,^(0))<|}, 

ns = {(i[o lfl> v r (r))<|}. 

Then limmf r _ >0 oP r (fli) > 1 ~ v/% since 5 < ft. The set {£ € M : 
(1[0,5 2 ])C) < e /3} is open in the weak topology. So by ([30]) . ([37]) . the Port- 
manteau theorem, and since 5 < 62, liminf r _ 5 . 00 P r (fi£) > l—rj/2. Similarly, 
by (j33|h (p8|) . the Portmanteau theorem, and since 5 < 63, liminf r _>. 00 P r (fig) 
1. Hence, 

liminf p r (n[ n fi£ n fiS) > 1 - ??. 

r—too 
Note that for all t > 0, at most one j € J(t) satisfies < Wj(t) < Wj by ([2]). 
So almost surely for all t > 0, SieJ(t) 1(0,5] ( w j(£)) 1(5, 00) (^j) < 1- Under 
fluid scaling, this implies 

sup (l m ,Z r (t)) < (l m ,Z r (0)) + (l m ,V r (T)) + -. 

tG[0,T] r 

Then almost surely on fi^ n fi?J ^ ^3> 

sup (l [0j(5] ,Z r (t)) <^ + ^ + -, 
te[0,T] 3 3 r 

which is bounded above by e for sufficiently large r. □ 

We are now ready to prove that the fluid scaled state descriptors under 
SJF and SRPT converge together. 

Theorem 2. For all T < 00, 

sup d[Z r (t), Z r Jt)\ A 0, as r ->■ 00. 
te[o,T] 
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Proof. Fix T < oo and e,r] > 0. By Lemma there exists (5 > such that 
the events 

ni = [ sup {(i m ,z r (t)) v (i m ,z;(t))) < f 1 

[te[o,T] ^J 

satisfy liminf r _ s . 00 P r (J-i) > 1 — r//2. By Lemma [61 the events 

n r 2 = \ SUp (Xl[ ,jrr (t) ), ££(£)) < — > 
[t_[0,T] 2 J 

satisfy liminfr-j-oo P r (il2) > 1 — '7/2, and so liminfr._j.oo P r (il^nil2) — 1 — '/■ 
Fix £ € [0, T] and let B C M+ be closed. By intersecting B with [0,5], 
(5,F r (t)), and [F r (t),oo), and by Markov's inequality, 

(IB, ^(O) < <1[0,5],^(0> + ^Xl[0,^(t))^ r W) + (l[F-(t)) )0 c)nB ^(t) 

We can replace Z r (t) by Z J p {i) almost surely in the last term by Corollary 
[21 and in the next to last term by Corollary [H since 

(Xl [0 ,^ (4 )),^(i)> = W r (t) - ( X l [P r {tloo) ,Z r (t)). 

Then use [F r (t)), oo) n B C B € to get, almost surely on __'J n D^, 



e e 

2 + 2 



<l B ,_*(t)><- + - + <l Be ,i*(t)). (39) 



The same argument yields, almost surely on £_^ D O^, 

(lB,Z;(t))<e+(l B *,Z r (t)). (40) 

Since t _ [0, T] and _? C R+ were arbitrary, (|39p and (|40p imply that, almost 
surely on HJ n ^j 

sup d[Z r (t),Z;(t)]<e. 

t€[0,T] 

D 

Proof of TheoremUl By (133H . -ZI(-) => _>*(•) as r — > oo. Theorem [T] thus fol- 
lows from Theorem [2] and the converging together lemma (see e.g. Theorem 
4.1 in [3J). □ 
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